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Abstract 

This paper reveals the relation between the canonical coset de- 
composition of unitary matrices and the corresponding decomposition 
via Householder reflections. These results can be used to parametrize 
unitary matrices via Householder reflections. 
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1 Introduction 



The parametrization of unitary matrices is important in many fields 
of physics such as quantum computation [TJ [2] , particle physics [3] and 
quantum iterferometry [4] . Several methods for the parametrization of 
SU(N) have been developed, including one that uses Euler angles [5] 
to find explicit expressions for the Haar measure. Another form was 
given by Reck et al, who explicitly constructed arbitrary unitary op- 
erators in terms of two-state pivots PQ. Dita proposed a parametriza- 
tion through diagonal unitary matrices interlaced with real orthogonal 
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matrices [6] and Rowe et al. developed another method studying the 
Wigner functions for SU(3) [3], which were also used in a formula- 
tion of coherent states for SU(N) in the work of Nemoto[7j- More 
recently, a technique resembling the canonical coset parametrization 
was presented by Jarlskog [H |9] . 

Matrix decomposition through Householder reflections |1U| was re- 
cently proposed as an efficient method for synthesizing unitary operators[2j 
111] . This technique can be implemented in certain quantum systems 
using N instead of N 2 steps as required in methods applying two- 
state pivots m following the logic of the Givens rotation |12] , The 
canonical coset parametrization |13t [H] of unitary matrices is a gen- 
eral method that may be used to calculate explicit expressions of the 
Haar [T3l El [15] and Bures measures [161 E] > which are known to be 
important in Bayesian quantum estimation [18J. The main purpose 
of this paper is to establish the connection between the Householder 
decomposition and the canonical coset parametrization of unitary op- 
erators, which could have practical significance in a number of physical 
applications. 



2 Householder Decomposition 

A unitary operator can be seen as a basis of orthonormal vectors and 
the Householder decomposition of unitary operators serves to align the 
original basis in terms of a new basis. The Householder decomposition 
consists of a sequence of transformations, each of which is a reflection 
with respect to a hyperplane defined by the orthonormal vector \n), 
such that 

R ln) =l-2\n)(n\. (1) 

As a proper reflection, Ri n \ satisfies det(R^) = — 1 and R\ n ^R\ n ^ = 1. 
Given a set of orthonormal basis elements je^), the unitary matrix 
U(N) may be decomposed as a product of N factors using Householder 
reflections that are designed to sequentially align the columns of U(N) 
along the orthonormal basis elements |e/t). 

In order to illustrate the Householder decomposition, let \Wi) be 
the first column of U(N) and 4>\ the phase of the topmost complex 
component of \W\). The first Householder reflection is constructed as 

R\ ui) = i_2— j—HmKuil, (2) 
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where \u±) = \Wi) + e 1 ^ 1 \ei) , such that i2i Ul \|Wi) = — |ei). For visual- 
ization, the case with two-dimensional real vectors is shown in Figure 
[TJ The most general definition of the vector \ui) requires \e\) to be 



\e 2 ) 



\ 1 











l«l> 



Figure 1: Graphical representation of a Householder reflection constructed 
for a real two-dimensional unit-length vector \Wi) reflected on the plane 
perpendicular to \ui). The net result is a reflection that positions \ Wi) along 
the direction — |ei). An additional reflection can be performed to position 
\Wi) along |ei). 



multiplied by the magnitude of the vector \W\), but as U(N) is uni- 
tary, then \Wi) has unit norm. Although, the opposite sign choice 
e wi _i. — e wi cou id be utilized, a positive sign is required in order 
to arrive at the canonical coset decomposition as shown in the proof 
of Theorem [T] below. Moreover, the positive sign is usually chosen 
to obtain better numerical stability (see [19], page 225). The above 
procedure is repeated in recursion, with IW2) as the second column of 
R\ U1 )U(N) (instead of U(n)) and \u 2 ) = \W 2 ) + e^ 2 |e 2 ). In this way, 
the resulting Householder decomposition becomes 

U(N) = R lui) R lu2) ...R luN _ l} U(l) N , (3) 

where U(1) N is the diagonal matrix with element e 1 ^ 1 at the l-th po- 
sition along the diagonal. 

The construct of the canonical coset decomposition from the House- 
holder decomposition is given by the following theorem 
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Theorem 1 The factors of the canonical coset decomposition 

U(N) = Wl U{N -V *W—U(1)°» 

{ ' U(N- 1)®U(1) U(N-2)® U{\) 17(1) <g>*7(l) 1 ; ' 

(4) 

can be expressed in terms of Householder reflections R\ Uk ) a $ 

C/(7V-1)®[/(1) " ^M^> ^ 
[7(iV - 2) " (6) 

17(1)®" = R lei) R le2) ...R leN _ l} U(l) N (7) 
where the corresponding Householder decomposition is 

U(N) = R lui) R lu2) ...R luN _ l} U(l) N , (8) 

such that 

\u 1 ) = \W 1 )+e i ^\e 1 ), (9) 

with \Wi) being the first column of U(N) and <f>\ the phase of the 
first component of \W\). The remaining vectors |n^) are calculated 
recursively. 

Proof. Let us define the normalized vector \n\) as 

The first column of U (N) can be decomposed into parallel and per- 
pendicular parts with respect to \e\) as 

\Wi) = \W n ) + \W 1± ) = /oe^lei) + \W 1± ), (11) 

where pe 1 ^ 1 is the polar representation of the first component of \W\). 
The normalized vector \n\) can be expressed as 

K) = J_ (\W 1± ) + (p + l)e^|ci>) • (12) 

This allows writing |ni) as 

\ n \) = \ n \\) + \ n ±) = 7|ei) + |n_L), (13) 
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with 7 = (^ +1 ) e and In i ) = ^ 1X ^ . From (1111 ) and (1121 ) we obtain 

VKW VMM 

1 = p 2 + (W 1± \W 1± ) (14) 
= (Wi±|P^i ± ) + (p+i) 2 , (15) 

which can be used to extract (ui\ui) = 2(1 + p) and write 7 as 



7 = yfi + i>«* (16) 
The product of the following Householder operators can be expanded 



as 



R\ n) R\ ei) = (l-2|n)(n|)(l-2|ei)( ei |) (17) 
= l-2|n ± )(n ± |-2 7 |ei)(n ± |+27> ± )(ei| 
+ 2(| 7 | 2 -l)| ei )( ei |. 

The identity matrix can be decomposed as 1 = lj_ + |ei)(ei| and 
introduced above to yield 

R\n)R\ ei ) = l±-2|n ± )(n ± |-2 7 |e 1 )(n ± |+27*|n±)(e 1 |+(2| 7 | 2 -l)|e 1 )(e 1 | 

(18) 

which can be decomposed into four parts expressed in block matrix 
form as 

R \n)%i) = (X*\n ± ) lj.-2K)(n ± |)- (19) 
Next, the following variable change is applied 

|X> = 27>j_>, (20) 
and the variable r is defined as the magnitude of \X) 



r = V(X\X). (21) 

From (|13p we see that 

r 2 =4| 7 | 2 (n ± |n ± )=4| 7 | 2 (l-| 7 | 2 ), (22) 
and with some algebra this leads to 

a/1 -r 2 = 2| 7 | 2 - 1. (23) 
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This construction is consistent only if 1 > 2 1 | 2 — 1 > 0, which is seen 
to be true by inspecting (fT6|) , Moreover, (fTB]) implies that 2 1 | 2 — 1 = p. 

The block matrix f)19|) can be written in the form of a canonical 
coset as follows 



R\n)R\ ei ) = |yx , L^.vwv,) ( 24 ) 



-(x\ 

\x) i L -i=^^\x){x\ t 

^ (25) 

, \X) y/l ± -\X)(X\) 

U(N) 



(26) 



C/(7V - 1) ® tf(l) ' 

which is the form given by Gilmore|131 114j. 

The remainder of the cosets can be found recursively with the 
insertion of reflections of the form R\ ek \, which can be constructed by 
replacing the k-th element of the identity matrix with —1. There is 
flexibility in the sequential choice of the reflections R\ ek ) because they 
commute with each other [R^,R^ e ^] = 0. Additionally, since R\ Uk ) 
is a block matrix, the following identity can be verified 

[R\ ej) ,R\ Uk )] = ioxk>j. (27) 
Thus, the theorem is proved. ■ 

Corollary 1 The factors of the reversed canonical coset decomposi- 
tion, 

y ' { ) U(l) ® 17(1) ' U(N - 2) ® U(l) U(N - 1) ® U(l) ' 

(28) 

can be expressed in terms of the reversed Householder reflections R\ Uk ) 
as 



U(N) 



U(N-1)®U(1) 

U(N-l) 
U{N - 2) i [/(!) 



Rle^m) (29) 

R \e 2 ) R \u 2 ) ( 30 ) 



17(1)®" = U(l) N R leN _ iy ..R ]e2) R ]ei) (31) 
where, the corresponding Householder decomposition is 

U(N) = U(l) N R luN _ l} ...R lu2) R lui) , (32) 
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such that 

(u 1 \ = (W 1 \ + e i ^(e 1 \, (33) 

with (Wi\ being the first row of W and (pi the phase of the first com- 
ponent of (Wi\. The remainder of the vectors {uk\ are calculated re- 
cursively. 

One important consequence of Theorem [1] is the opportunity it af- 
fords to parametrize the Householder reflections. This can be accom- 
plished with the use of ([20)) . as described in the following corollary. 



Corollary 2 The normal Householder vector \n) can be parametrized 
in terms of the canonical coset vector \X) as 

\n)= 1 \e 1 ) + ^\X), (34) 

where 

l7l= ^P> +1 . (35, 

with the phase of 7 extracted from the phase of the first diagonal ele- 
ment ofU{l) N . 

The parametrization of unitary matrices is essential in order to 
calculate the Haar metric and Haar measure [15J of unitary matrices, 
and for this purpose using Householder reflections is more efficient 
than other alternatives such as the parametrization in terms of Eu- 
ler angles [5]. In the same spirit, the parametrization of the House- 
holder normal vectors \n) (134p can be exploited to generate random 
unitary matrices with the Haar measure as suggested by Ivanov |20j . 
This can be accomplished by homogeneously distributing \X) within 
Euclidean balls [21J while uniformly distributing the phase of 7 in 
the range arg{^) G [— ir,ir}. For example, the generation of random 
unitary matrices U(3) with the Haar measure requires the following 
parametrization 

17(3) = (1 - 2|m)(ni|)(l - 2\n 2 ){n 2 \)diag{e i 't>\e i * 2 ,e^ 3 ), (36) 

where the phases (j>k are drawn from a uniform distribution over the 
range [— vr,7r] and where \n\) and \ri2) are parametrized within balls 
B A and B 2 , respectively. 
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As an illustrative example for the calculation of the canonical coset 
decomposition, consider 



/ i i 

/ ./9 ./O 



s/2 \/2 

i i 

"2 2 

1 1 

"2 2 







1 



(37) 



The Householder decomposition is Uq = |i?)| Ul ^|i?)| M2 ^(l) 3 , where 



l" R )|«2) 



W(1) S 



1 

2 

(2+t/2) 
4 



4+2^ 
2+2^ 
4 



\ 2 4+2^ 

A o o 

_ J_ i(l+%/2) 



n «(l+>/2) 
2+^2 




2+V2 
1 



(38) 



(39) 



(40) 



This leads to the corresponding canonical coset decomposition 

1 \ 

9 \ 



U(3) 



17(2) ® 17(1) 

g(g) 
17(1)® C/(l) 




4+2^ 
2+2y^ 
4 



\ 





z(1+a/2) 
2+^/2 

73 / 



(41) 



(42) 



(43) 



The reversed canonical coset decomposition (j28|) can be obtained 



<s 



using the reversed Householder reflections with the following results 



17(3) 



17(2) ® U(l) 
U(1)®U(1) 




(44) 

(45) 
(46) 



3 Conclusions 

We formally proved the connection between the Householder and canon- 
ical coset decompositions of unitary matrices. This result allows for 
performing the canonical coset decomposition of unitary matrices very 
efficiently. Furthermore, the Householder decomposition is now fully 
parametrized in terms of the canonical coset vectors \X) (I34p . 
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Appendix 



The following exponential involving the complex vector B can be writ- 
ten in terms of the Cartesian coordinates x 3 as 



cxp 



B 

-St 



-7ft [1-X^X] 1 /' 2 



(47) 



such that 



/ x 1 + ix 2 



B^B 



X + IX 



2N-3 + ix 2N-2j 



(48) 
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The Cartesian coordinates range inside an even ball B 2k , where the 
radial coordinate is r 2 = X<X. This exponential is important because 
it provides a parametrization of the coset [/(jv-i)2t/(i) as a iV x A?" 
matrix, which can be used to parametrize the unitary operator Q, as 

^ G U(N) U(N) U(N-l) U(2) 



U(1)® N U(N-1)®U(1) U(N-2)<8>U(l)'"U(l)<8>U(iy 

(49) 

In the current paper we are interested in the representation of U(3), 
which are generated from two cosets. The left coset is 



U(2) 







U(1)®U(1) 



a/ 1 - ( x1 ) 2 - ( x2 ) 2 ~ xl + ix2 

K o ' x 1 + ix 2 V 1 - ( xl Y - (^ 2 ) 2 / 

(50) 

with the variables inside a disk (x 1 ) 2 + (x 2 ) 2 < 1. The right coset is 



U(3) [V^F ~x Fj + ix 6 -x 3 + ix 4 ^ 

w w \x 3 + zx 4 y 2 * 3 v 33 



with 



y 22 = (x 3 ) 2 + (x 4 ) 2 + y^^((x 5 ) 2 + (x 6 ) 2 ) (52) 



V 23 = V ' 1 £ 1 (x 3 -ix 4 )(x 5 + ix 6 ) (53) 



^33 = v/ y^ (x 3 + x 4 ) + l(x 5 + x 6 ) (54) 

^2 = (x 3 )2 + (x 4 )2 + (x 5 } 2 + (x 6 )2 (55) 

and the variables in the following range 

(x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 + (x 6 ) 2 < 1 (56) 
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